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Abstract 

Hard-photon (E^ > 30 MeV) emission from second-chance nucleon-nucleon Bremsstrahlung 
collisions in intermediate energy heavy-ion reactions is studied employing a realistic 
thermal model. Photon spectra and yields measured in several nucleus-nucleus re- 
actions are consistent with an emission from hot nuclear systems with temperatures 
T ~ 4 - 7 MeV. The corresponding caloric curve in the region of excitation energies 
e* ~ 3^ - 8^ MeV shows lower values of T than those expected for a Fermi fluid. 
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1 Introduction 



The determination of the thermodynamical properties such as temperature, 
density, and excitation energy of the hot nuclear systems produced in nucleus- 
nucleus reactions is one of the main goals of heavy-ion (HI) physics. At mod- 
erate excitation energies, e* ?a 3^4 - 15^ MeV, the experimental derivation of 
these observables as well as their correlation is a prerequisite for a quantitative 
investigation of the nuclear equation of state e = e(p,T) in connection with a 
possible liquid-gas phase transition [1-4]. To date the most unambiguous evi- 
dence for such a phase transition in HI collisions is given by the "caloric curve" 
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[1] which relates the thermal energy of an excited nucleus to its temperature, 
e* = e*(T). For very low e*, the caloric curve of nuclei is found to follow very 
closely [5] the ~ Fermi law given by the well-known Bethe formula for the 
nuclear density of states [6] . Such a trend is maintained at moderately low nu- 
clear excitation energies (e* <3^ MeV) [7]. In the region e* 3^ - 8^ MeV, 
however, the curve flattens [1] suggesting a phase transition (the width of the 
"plateau" indicating the latent heat related to the phase change). Although 
such thermodynamical behaviour has also been observed in other microscopic 
systems such as metallic [8] and hydrogen [9] clusters, the empirical determina- 
tion (and the theoretical interpretation) of the nuclear caloric curve has been 
a much debated issue the last years [1-4,10,11]. As a matter of fact, the three 
experimental methods employed so far to measure the temperatures attained 
in a reaction do not yield fully equivalent caloric curves. The existing nuclear 
"thermometers" are based on: (i) the slopes of the kinetic energy spectra of 
light particles (n, p, a) [12], (ii) the double ratios of neighbour isotopes [13], 
and (iii) the relative populations of excited states [14]. Having in hand an 
alternative thermometer based on a clean and weakly interacting probe would 
be extremely useful when searching for signals of the nuclear liquid-gas phase 
transition. Electromagnetic probes, viz. photons and dileptons, due to their 
weak final state interaction with the surrounding medium, have long been 
recognized as the most direct probes of the space-time evolution of the collid- 
ing nucleons [15]. In HI reactions at intermediate bombarding energies, two 
concurrent processes arc known to be responsible for hard-photon emission 
[16]: (i) first-chance (off-equilibrium) collisions between projectile and target 
nucleons in the first stages of the reaction, and (ii) subsequent NN coUisions 
in the produced "nuclear fireball" zone. The experimental existence of a sec- 
ond thermal component emitted from the produced hot nuclear sources and 
accounting for up to 30% of the total photon yield above E^=30 MeV, has 
been unambiguously observed recently [17-20]. In this Letter we propose to 
exploit such radiation as a novel thermometer of hot nuclear matter using a 
realistic thermal model which reproduces satisfactorily the observed photon 
spectral shapes and yields. Using such a thermometer we then construct the 
caloric curve in the region of the expected liquid- gas phase change. 



2 Thermal model 



Photons produced in HI reactions escape freely the interaction region imme- 
diately after their production. Thus, even when emitted from an equilibrated 

source they do not have a blackbody spectrum at the source temperature. 
However, the inverse slope parameter of their spectrum Eq and the tempera- 
ture T of the nuclear medium are strongly correlated. In this work we employ 
the thermal model of Neuhauser and Koonin [21] (henceforth NK) to quan- 
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titatively relate and T. According to this model, the differential rate of 
photons emitted in incoherent NN^ processes within a hot nuclear fragment 
is described by the expression: 



d^xdtdE, (2^(2;0^^^P^'^^^P^'^^^^'dK;' 



where pi^2i and Pi2i are the initial momenta and relative- velocity of the collid- 
ing nucleons, /(p) their (single-particle) momentum distribution, and da^/dE^ 
the angle-integrated Pauli-blocked NN Bremsstrahlung cross-section. Approx- 
imating the emitting region as nuclear matter in thermal equilibrium with 
local temperature T and density p, the momentum distribution can be simply 
parametrized by a hot Fermi-Dirac distribution {h — c — — 1)' 



/(P) 



1 -|- exp 



/T 



(2) 



normalized to 4 / dp/(p)/(27r)^ = p. For low temperatures (T <g ep) the 
chemical potential // of a system of nucleons can be written as a function of 
the Fermi energy ei?(p) [where ep — \Jpf + mj^—mN andpi?(p) = (37r^p/2)^/^] 
and T [22]: 



T 



^f{p) 



(3) 



The in-medium Bremsstrahlung cross-section da^/dE^ of Eq. (1) is approxi- 
mately [21] 

da^ da^, r dfl^dflf 



dE^ dE^ 



/ll-/(Pi/)l[l-/(P2,)1^^5i/, (4, 



where da^/dE^ is the elementary A^A^ Bremsstrahlung cross-section in free 
space, [1 — /] the usual Pauli-blocking factors, and dQ^j the solid angle of 
the outgoing gamma and nucleons. At the considered energies, the reaction 
pp — > PP7 is suppressed by a factor ~20 with respect to pn — > pn'j [23], 
neutron-neutron Bremsstrahlung is vanishingly small, and thus one needs only 
to consider the pn^ process. The isospin- averaged cross-section da^/dE^ is 
then one half of da-pn^/ dE^. We employ here the parametrization of Schafer 
et al. for dcjpn-y/ dE^ derived within a proper relativistic and gauge-invariant 
meson-exchange effective model for the A^A^ interaction [23], which reproduces 
well the available data on pn Bremsstrahlung [24] . 
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2.1 Nuclear temperatures 



The emission rates given by Eq. (1) can be thus calculated for a nuclear system 
at temperature T and density p. The Bremsstrahlung rates are shown in Fig. 
1 for a source at T = 4, 6 MeV and p/po = 0.5, 1.0 (po= 0.16 fm~'^). The 
resulting hard-photon distributions above 30 MeV can be well approximated 
by a Boltzmann exponential with slope Eq in agreement with the experimental 
data (Fig. 2). The integrated yields scale approximately with T^ '' and p. The 
temperature T of the emitting source and the photon slope parameter £^q, 
extracted from an exponential fit of the model spectra above = 30 MeV, 
are found to be well described by the relation: 

T(MeV) = (0.78 ± 0.02) • ^o(MeV), (5) 
in the range T ^ 3 - 10 MeV and p ^ (0.3 - 1.2)po. 

Applying Eq. (5) the nuclear temperature attained in a HI reaction can be 
determined by measuring the slope of its thermal hard-photon spectrum. The 
temperatures obtained for all the systems studied by the TAPS collabora- 
tion where a thermal hard-photon component has been measured [17-20] lie 
in the range T ^ A - 7 MeV (Table 1) in agreement with the typical values 
found in intermediate-energy nucleus-nucleus collisions [2,3]. The highest T 
(Eq) values for a given incident energy eiab correspond to the most symmet- 
ric systems, which have the highest energy deposition in the center-of-mass, 
eAA = Ared^iab/ Atot, i-G., the largest attainable excitation energies. In compar- 
ison with the nuclear thermometers used so far [12-14], the thermal-photon 
thermometer presents several advantages: (i) minimal preequilibrium contam- 
ination due to the large difference, by a factor two to three, between the direct 
(first-chance) and thermal hard-photon slopes [18]; (ii) absence of final-state 
distortions like side feeding, rescattering, and reaccelaration by the Coulomb 
field and/or collective motion; (iii) measurement of the T of the system right 
after equilibration, just before its break-up [19], when the maximum thermal 
energy of the equilibrated residue is achieved; and (iv) intrinsic selection of 
(semi)central dissipative reactions with large number of A^A^ collisions [19]. 
The use of thermal photon slopes as a nuclear thermometer is however re- 
stricted to intermediate-energy reactions with eiab ~ 30^-90/1 MeV. The lower 
limit in eiat results from the experimental difficulty to resolve accurately the 
thermal component in a double-exponential fit of the hard-photon spectrum 
due to: (i) the very small cross-section in the direct high-energy regioriQ and 
(ii) the increasing role at E^ < 20 MeV of statistical 7 from decays of Giant 
Dipole Resonances and bound states. Above eiab ~ 90^ MeV, hard-photon 

^ The emission probability of a photon of E^y > 30 MeV per pn colhsion is a steeply 
increasing function of the incident energy ei^b [18]: ~ exp[— (1/eiab)]- 
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spectra are well described by a single "direct" exponential [29,30] and also the 
background of photons from the decay of the produced vr^'s becomes important 
[30]. 

2.2 Photon multiplicities 

The absolute thermal photon yield per nuclear reaction, = ct^/o-aa, pre- 
dicted by the NK model can be obtained integrating Eq. (1) over the relevant 
space-time history of the equilibrated system produced in a nucleus-nucleus 
reaction: 

Ej=30 MeV ' 

In the most general case T = T{x,t) and p = p{x,t), and the integration 
in Eq. (6) requires a consistent modeling of the space-time evolution of the 
reaction as done, e.g., in hydrodynamical approaches. Since we are just inter- 
ested in verifying that the thermal model provides also a correct estimation 
of the measured thermal photon yields, we will significantly simplify the cal- 
culation of M^^ making a few plausible assumptions. First, one may neglect 
any T and p gradients within the nuclear source and approximate the inte- 
gral over space by the volume of the participant nucleons in the reaction, 
V = Atot/p, where Atot = Ai + A2 is the sum of target and projectile nucle- 
ons[^. Secondly, one can simply replace the integral over time by a constant 
interval equal to the lifetime of the radiating source. At ^ 100 fm/c, consis- 
tent with the usual estimates of the lifetime of equilibrated excited nuclei [3]. 
With these rough approximations, partial integration of Eq. (1) above E^= 
30 MeV permits to compute M^^ straightforwardly. Since the emission rates 
given by Eq. (1) scale as ~ p and since V oc p~^, the thermal 7 multiplicities 
are rather insensitive to the density of the source in this very simple estima- 
tion. Most experimental multiplicities (Table 1) are very well reproduced by 
the integrated NK photon rates and our simplified ansatz of the space-time 
history of the radiating sourcej^. The photon yields of the two systems with 

^ The equilibrated nuclear object produced in (semi) central collisions at these en- 
ergies is likely to include many more constituents than just the original partici- 
pants as the excitation is shared with other nucleons. Although it is improbable 
that all nucleons are available to participate in the nuclear Bremsstrahlung process 
(e.g. 10%-20% of them may just be gone in preequilibrium emission), we take here 
Atot = ^1 + ^2 as a reliable measure of the (relative) size of the system in the 
different reactions. 

^ The reported errors on M^^ in Table 1 include the uncertainties in T propagated 
from Eq and are quite large due to the strong power dependence of the photon 
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largest excitation energies (^^Kr+^^Ni and ^^Ar+^^Ni at 60^ MeV), however, 
are overpredicted by the model (though still within the combined M^^ and 
M™p systematic uncertainties). A much better agreement would be reached 
in these two cases if the lifetime of the equilibrated sources was reduced by a 
factor ~2 compared to lower energies (i.e. if Ar ~ 50 fm/c). This result is in 
quantitative agreement with the observed fast break-up of nuclear systems for 
increasingly larger values of e* interpreted in terms of a low-density (spinodal) 
fragmentation [31]. 



3 Caloric curve 

Having determined the values of T from through Eq. (5) for our different 
reactions, a caloric curve e*(T), can be constructed correlating T with the 
thermal excitation energies e* attained in each reaction. To construct such a 
curve, we use the published values of e* measured in semi-central and central 
reactions with equivalent systems (Table 1, most right column). In the case 
of ^^Kr-|-^^Ni at 60^4, in the absence of an experimental measurement of e*, 
an upper value can be obtained from the total (Coulomb-corrected) center- 
of-mass energy, caa, subtracted of other energy contributions^. Although it 
is fair to acknowledge that our globally indirect determination of e* is po- 
tentially subject to significant systematic uncertainties, such uncertainties are 
reasonably contained in the already intrinsically large experimental errors of 
the excitation energy measurements [2] and, furthermore, the qualitative con- 
clusions that we obtain from the analysis of the photon caloric curve would 
only be modified by much larger shifts along the e* axis. The e* - T correla- 
tion for the 7 reactions considered here is shown as solid dots in Fig. 3. Our 
caloric curve falls somewhat above the ALADIN [1] and EOS [32] resultsQ 
obtained with two different isotopic ratios, but still clearly below the region 
expected for a pure degenerate Fermi fluid (Fig. 3, dark band). The fact that 
the thermal photon temperatures seem to be slightly larger than the isotopic 
ones is not surprising since "chemical" temperatures probe a later lower den- 
sity (cooler) stage when isotopes are produced, whereas the most significant 
part of thermal photon emission takes place somewhat earlier. Indeed in this 



emission rates on T. 

^ e* (fiAA + Q) — Ccoii — fiprocq — Erot, wheiei (i) the reaction Q-value and e^ot are of the 
order MeV and have opposite signs, and (ii) the radial flow energy, ecoii, and 
the pre-equilibrium component, eprccqi represent r^3A MeV for a symmetric system 
at these energies [1~3]. 

^ Note however that both, ALADIN and EOS, collaborations have updated slightly 
upwards their original curves (shown in Fig. 3) to account for several corrections 
in T (sequential feedings) and e* (collective radial flow energy), see e.g. [11] for a 
compilation of the most recent results. 
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scenario, the radiation of thermal Bremsstrahlung photons occurs during the 
path followed by the excited nuclear system between the hot "liquid" phase 
(one single big nuclear fragment at the end of the pure pre-equilibrium stage) 
and the mixed phase of gaseous particles and small droplets produced later in 
the fragmentation process [19]. Our caloric curve disagrees, on the other hand, 
with the one obtained by the INDRA collaboration using kinetic temperatures 
[10], which follows closely the e* ~ trend. The kinetic T's are systemat- 
ically higher by about 1-2 MeV most likely due to the fact that the slope 
parameters of p and a kinetic energy spectra are influenced by preequilibrium 
and mid-rapidity emissions [33], and/or nucleon Fermi motion effects. The de- 
parture of the photon caloric curve from the e* = aT^ law (with level density 
parameters a = A/8 — A/13 MeV~^) characteristic of the nuclear ground state 
(liquid phase) , and the slow increase of T with e* in the region of intermediate 
excitation energies are in qualitative agreement with the expected behaviour 
of a liquid-gas phase transition occurring in excited atomic nuclei. 



4 Conclusions 

A thermal bremsstrahlung model has been employed to extract the thermody- 
namical properties of the nuclear systems produced in heavy-ion collisions at 
intermediate-energies. Such a model predicts exponential hard-photon spectra 
above — 30 MeV in agreement with the experimental data. The thermal 
slopes are linearly correlated with the temperature of the emitting system. 
This defines a new thermometer of hot nuclear matter which, at variance with 
the usual hadronic-based methods, is free of significant final-state distortions. 
The hot nuclear residues prepared in different reactions with excitation ener- 
gies e* ~ 3^ - 8^4 MeV, have temperatures in the range T f» 4 - 7 MeV yielding 
a slowly increasing caloric curve which hes below the expected behaviour for 
a pure Fermi liquid. 
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Fig. 1. Thermal Bremsstrahlung emission rates, Eq. (1), given by the NK model 
for a nuclear system in equilibrium at temperatures T = 4, 6 MeV and densities 
p/po=0.5, 1.0. 
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Fig. 2. Thermal hard-photon spectrum measured in the reaction ^^Ar+^^^Au at 60^ 
MeV [19] compared (dashed band) to the prediction of the NK model for a source 
at p = po and r=5.3 ± 0.5 MeV (see text and Table 1). 
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Fig. 3. Caloric curve constructed with the photon slope thermometer (dots) com- 
pared to ALADIN (squares) and EOS (crosses) curves (isotopic T's), and to INDRA 
(rhombbi) curve (kinetic T's). The dark area corresponds to the Fermi relation 
= aT^ with a = A/S- A/13 MeV-^ 
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Table 1 

Heavy-ion reactions studied by the TAPS collaboration where a thermal 
bremsstrahlung component has been identified. For each reaction we report: (i) 
the nuclear temperatures T extracted through Eq. (5) from the measured thermal 
slopes Eq-, (ii) the total number of nucleons Atot, the photon multiplicity M^^ pre- 
dicted by the NK model (for a source with volume V = Atot/ Pi and lifetime At=100 
fm/c) compared to the experimental value M^""; and (iii) the excitation energies e* 
measured in other experiments. 



System 


e,,b (^MeV) 


(MeV) 


T (MeV) 


^tot 


(10-4) 




e* (^MeV) 


208pb+197Au 


30 


5.5 ± 0.6 [17] 


4.3 ± 0.5 


405 


0.7 ± 0.7 


0.6 ± 0.2 [17] 


3.7 ± 0.7 [25] 


36Ar+197Au 


60 


6.8 ± 0.6 [19] 


5.3 ± 0.5 


233 


1.9 ± 1.3 


1.6 ± 0.2 [19] 


5.1 ± 1.0 [26] 


i8iTa+i97Au 


40 


6.9 ± 0.6 [17] 


5.4 ± 0.5 


378 


3.4 ± 1.9 


3.2 ± 1.0 [17] 


7.0 ± 1.5 [27] 


36Ar+i07Ag 


60 


7.0 ± 1.0 [19] 


5.5 ± 0.8 


143 


1.5 ± 1.5 


1.2 ± 0.2 [19] 


6.3 ± 1.2 [26] 


I29xe+ii2sn 


50 


7.0± 0.6 [20] 


5.5 ± 0.6 


241 


3.1 ± 2.2 


2.6 ± 0.4 [20] 


7.0 ± 1.5 [28] 




60 


8.5 ± 0.8 [17] 


6.6 ± 0.6 


144 


4.4 ± 2.7 


2.0 ± 0.4 [17] 


<10. (see text) 


36Ar+58Ni 


60 


8.8 ± 1.0 [19] 


6.9 ± 0.8 


94 


4.1 ± 2.9 


1.1 ± 0.2 [19] 


8.0 ± 1.5 [26] 
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